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Abstract. A linear system of differential equations describing a joint motion of 
thermoelastic porous body and thermofluid occupying porous space is considered. 
Although the problem is linear, it is very hard to tackle due to the fact that its 
main differential equations involve non-smooth oscillatory coefficients, both big 
and small, under the differentiation operators. The rigorous justification, under 
various conditions imposed on physical parameters, is fulfilled for homogenization 
procedures as the dimensionless size of the pores tends to zero, while the porous 
body is geometrically periodic. As the results, we derive Biot's like system of 
equations of thermo-poroelasticity, system of equations of thermo-viscoelasticity, 
or decoupled system consisting of non-isotropic Lame's equations for thermoelastic 
solid and Darcy's system of filtration for thermofluid, depending on ratios between 
physical parameters. The proofs are based on Nguetseng's two-scale convergence 
method of homogenization in periodic structures. 

Key words: Biot's equations, Stokes equations, Lame's equations, two-scale con- 
vergence, homogenization of periodic structures, thermo-poroelasticity, thermo- 
visco-elasticity. 

Introduction 

In this article a problem of modelling of small perturbations in thermoe- 
lastic deformable solid, perforated by a system of channels (pores) filled with 
thermofluid, is considered. The solid component of such a medium has a 
name of skeleton, and the domain, which is filled with a fluid, is named a 
porous space. The exact mathematical model of such of medium consists of 
the classical equations of momentum, energy and mass balance, which are 
stated in Euler variables, of the equations determining stress fields and ther- 
modynamics law in both solid and liquid phases, and of an endowing relation 
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determining behavior of the interface between liquid and solid components. 
The latter relation expresses the fact that the interface is a material surface, 
which amounts to the condition that it consists of the same material parti- 
cles all the time. Clearly the above stated original model is a model with 
an unknown (free) boundary. The more precise formulation of the nonlinear 
problem is not in focus of our present work. Instead, we aim to study the 
problem, linearized at the rest state. In continuum mechanics the methods 
of linearization are developed rather deeply. The so obtained linear model 
has been studied in [11]. Further we refer to this model as to model (NA). 
In this model the characteristic function of the porous space x is a known 
function for t > 0. It is assumed that this function coincides with the char- 
acteristic function of the porous space Xi given at the initial moment. 
In dimensionless variables (without primes) 



the differential equations of model (NA) in the domain fteK 3 for the 
dimensionless displacement vector w of the continuum medium and the di- 
mensionless temperature 9, have the form: 



x' = Lx, t' = rt, w' = 



Lw, 9' = ®*—9 



TV* 



d 2 w 



div x ¥ + pF, 



(0.1) 



3t 2 




(0.2) 



(0.3) 



(0.4) 




p + xa P div x w = 0, 



(0.5) 
(0.6) 



Here and further we use notations 



D(x, u) 



(1/2) (V x u + {V x u) T ) , 



P = XPf + (1 
«x = + (1 



X)Ps, c p = xcpf + (1 - x)c ps , 
X)a„ s , a e = xa ef + (1 - x)®e s 
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The dimensionless vector of distributed mass forces and dimensionless 
mean densities of the solid and the fluid at rest (all without primes), respec- 
tively, are given by the formulas 



F' = gF, p' s = pop,, p' f = 
and the dimensionless ratios by the formulas 

2/1 



a T - 

OLv ~- 

ocx 

otes 

c Pf 



9 r 

(v ~ f P) 

vlrgpoL ' 
2\ 



a. 



a. 



rgpoL 

3 



a 



gpoL 

7sP^* 



a. 



v*rgp 
Cfmp'f'&l 

Povl 



aof 



ps 



gpoL 

XfV* 

vlrgpoV 
( 2c f P + p'f c f PP )(p'f 

gpoL 

v*rgp 
Pavi 



I \2 



(0.7) 



where all dimensionless constants on the left-hand sides of relations are pos- 
itive and all constants on the right-hand sides of relations are described in 

m 



We endow model (NA) with initial data 

dw , 



W\ t=0 = Wq, 



dt 



\t=0 



«0> 



U=o 



and homogeneous boundary conditions 

0, 9 = 0, xeS = dtt, t>0. 



w 



(0.8) 



(0.9) 



From the purely mathematical point of view, the corresponding initial- 
boundary value problem for model (NA) is well-posed in the sense that it 
has a unique solution belonging to a suitable functional space on any finite 
temporal interval (see [H]). However, in view of possible applications, for 
example, for developing numerical codes, this model is ineffective due to its 
sophistication even if a modern supercomputer is available. The differential 
equations of the model involve frequently oscillating non-smooth coefficients, 
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which have structures of linear combinations of the function % These coeffi- 
cients undergo differentiation with respect to x and besides may be very big 
or very small quantities as compared to the main small parameter e. In the 
model under consideration we define the dimensionless size of the pores e as 
the characteristic size of pores I divided by the characteristic size L of the 
entire porous body: 

I 

£ = T 

Therefore a question of finding an effective approximate models is vital. Since 
the model involves the small parameter e, the most natural approach to this 
problem is to derive models that would describe limiting regimes arising as 
e tends to zero. Such an approximation significantly simplifies the original 
problem and at the same time preserves all of its main features. But even 
this approach is too hard to work out, and some additional simplifying as- 
sumptions are necessary. In terms of geometrical properties of the medium, 
the most appropriate is to simplify the problem postulating that the porous 
structure is periodic. Further by model (NB) £ we will call model NA 
supplemented by this periodicity condition. Thus, our main goal now is a 
derivation of all possible homogenized equations in the model (NB) e . 

Remark 0.1. The particular case of the model (NA) when the system (10. ip - 
(10.61) is decoupled (agf = ag s = 0) we call as model (A). Correspondingly, 
the same particular case of the model (NB) £ we call as model (B) £ . 

The first research with the aim of finding limiting regimes in the case 
when the skeleton was assumed to be an absolutely rigid isothermal body was 
carried out by E. Sanchez-Palencia and L. Tartar. E. Sanchez-Palencia [T4"t 
Sec. 7.2] formally obtained Darcy's law of filtration using the method of two- 
scale asymptotic expansions, and L. Tartar [T4"l Appendix] mathematically 
rigorously justified the homogenization procedure. Using the same method of 
two-scale expansions J. Keller and R. Burridge |3j derived formally the system 
of Biot's equations [I] from the model (B) £ in the case when the parameter 
a M was of order e 2 , and the rest of the coefficients were fixed independent of 
e. It is well-known that the various modifications of Biot's model are bases 
of seismic acoustics problems up-to-date. This fact emphasizes importance 
of comprehensive study of model (NA) and model (NB) £ one more time. J. 
Keller and R. Burridge also considered model (B) e under assumption that 
all the physical parameters were fixed independent of e, and formally derived 
as the result a system of equations of viscoelasticity. 
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Under the same assumptions as in the article [3], the rigorous justification 
of Biot's model was given by G. Nguetseng [T3] and later by A. Mikelic, R. 
P. Gilbert, Th. Clopeaut, and J. L. Ferrin in O El E] - Also A. Mikelic et al 
derived a system of equations of viscoelasticity, when all the physical param- 
eters were fixed independent of e. The most general case of the model B e 
has been studied in [10]. In these works, Nguetseng's two-scale convergence 
method [HI EE2] was the main method of investigation of the model (B) £ . 

In the present work by means of the same method we investigate all 
possible limiting regimes in the model (NB) £ . This method in rather simple 
form discovers the structure of the weak limit of a sequence {z 6 } as e \ 0, 
where z e = u e v e and sequences {u 6 } and {v £ } converge as e \ merely 
weakly, but at the same time function u e has the special structure u e (x) = 
u(x/e) with u{y) being periodic in y. 

Moreover, this method allows to establish asymptotic expansions of a 
solution of the model (NB) e in the form 

w e (x, t) = e 13 (w (x, t) + ewx (x, t, x/e) + o(e)) , (0.10) 

where Wq(x, t) is a solution of the homogenized (limiting) problem, Wi(x, t, y) 
is a solution of some initial-boundary value problem posed on the generic pe- 
riodic cell of the porous space, and exponent (3 is defined by dimensionless 
parameters of the model. In some situations expansion (10.101) has more com- 
plicated form like 

w e (x, t) = e^ f (wl(x, t) + ew[ (x, t, x/e) + o(ej\ , 

in the liquid component and 

w £ (x, t) = e' 3 " (w s (x, t) + ew{ (x, t, x/e) + o(e)) , 

in the rigid component. For example, if a M = e 2 and ot\ = e" 1 , then j3f — — 1, 
wl = and /3 S = 1. 

Suppose that all dimensionless parameters depend on the small parameter 
e and there exist limits (finite or infinite) 

lim aJs) — Ho, limaA(e) = Ao, lima T (e) = r . 

e\0 e\0 e\0 

We restrict our consideration by the cases when when r < oo and either of 



5 



the following situations has place: 

(I) fjto = 0, < A < oo; 

(II) < fi < oo, A = oo; 

(III) < /i , A < oo. 

If To = oo then, re-normalizing the displacement vector and temperature by 
setting 

w — > a T w, 9 — > a T 8, 

we reduce the problem to one of the cases (I)-(III). 

In the present paper we show that for the first case the limiting regime 
is a two-velocity and two-temperature continuum media, which is described 
by different types of Biot's like system of equations of poro-elasticity coupled 
with corresponding heat equations (equations of thermo-poro-elasticity) , or 
one-velocity and one-temperature continuum media, which is described by 
non-isotropic Lam'e's system of equations coupled with corresponding heat 
equation (equations of thermo-elasticity) (theorem I2.2H . For the second case 
the limiting regime is a two-temperature continuum media, which is described 
as a first approximation by different types of Darcy's like system of equations 
of filtration for the velocity of the liquid component, coupled with correspond- 
ing heat equations (for this first approximation the solid component behaves 
yourself as absolutely rigid body) and as a second approximation - non- 
isotropic Lame's system of elasticity for the re-normalized displacements of 
the solid component, or two- velocity and two-temperature continuum media, 
which is described by Biot's like system of equations of thermo-poro-elasticity 
for the re- normalized displacements (theorem I2.3p . Finally, for the last case 
the limiting regime is a one- velocity and one-temperature continuum media, 
which is described by system of non-local equations of thermo-visco-elasticity 
or non-local equations of thermo-elasticityftheorem I2.4p . 

§1. Model (NB) £ 

1.1. Geometry of porous space. In the model (NB) £ the porous 
space has geometrically periodic structure. Its formal description is as follows 
0(2]. 

Firstly a geometric structure inside a pattern unit cell Y = (0, l) 3 is 
defined. Let Y s be a 'solid part' of the cell Y. The 'liquid part' Yf is its open 
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complement. Set Y]f := Y s + k, k G Z 3 , the translation of Y s on an integer- 
valued vector k. Union of such translations along all k, E s : = eZ3 Y^ is 
the 1-periodic repetition of Y s all over M. 3 . Let Efbe the open complement of 
E s in IR 3 . The following assumptions on geometry of Yf and Ef are accepted. 

(i) Y s is an open connected set of strictly positive measure with a Lipschitz 
boundary, and Yf also has strictly positive measure on Y. 

(ii) Ef and E s are open sets with C 0,1 -smooth boundaries. The set Ef is 
locally situated on one side of the boundary dEf, and the set E s is locally 
situated on one side of the boundary dE s and connected. 

Domains Q £ s and f2j are intersections of the domain Q with the sets M 3 
and M 3 , where the sets M 3 and M 3 are periodic domains in K 3 with generic 
cells eY s and eYf of the diameter e, respectively. 

Union Y s U F y is the closed cube Y = {y Gl 3 , < y t < 1, i = 1, 2, 3}, 
and the interface T £ = 9^2^ n9r2j is the ^-periodic repetition of the boundary 
£7 = sdYf fl dY s all over cu. 

Further by x — X £ we wiU denote the characteristic function of the porous 
space. 

For simplicity we accept the following constraint on the domain Q and 
the parameter e. 

Assumption 1.1. Domain Q is cube, Q := (0, l) 3 , and quantity 1/e is 
integer, so that Q always contains an integer number of elementary cells Yf. 

Under this assumption, we have 

X(x) = X £ (x) = X (x/e), (1.1) 

c p = cf p (x) = x £ (x)c pf + (1 - x £ (x))c ps , 

P = P £ ( x ) = X £ (x)p f + (1 - X £ (x))p s , 

a„ = a%(x) = x £ (x)a„ f + (1 - x £ {x))a„ s , 

a e = a e (x) = x £ (x)a e f + (1 - x £ (x))a es , 

where x{v) is the characteristic function of Yf in Y . 

We say that a porous space is disconnected (isolated pores) if 

7 n = 0. 

1.2. Generalized solutions in the model (NB) £ . Equations (10. ip - 
(10.61) are understood in the sense of distributions theory. They involve the 
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equations (10. ip ~ (10. 6p in a usual sense in the domains Q'j and Q £ s and the 
boundary conditions 

[0] = 0, [to] = 0, x E T £ , t > 0, (1.2) 
[P] = o, [a%V x 6] =0, x E T £ , t > (1.3) 

on the interface T £ , where 

[(p](x ) = (p{ s )(x ) - (P(J)(X ), 
(P( s )(x ) = lim (p(x), if(f)(x ) = lim ip(x). 



x — ► aso x — • xq 

x G n't x g Q j 



There are various equivalent in the sense of distributions forms of rep- 
resentation of equations (I0.1I) - (I0.2I) and boundary conditions (11.2D — (11.3D . In 
what follows, it is convenient to write them in the form of the integral equal- 
ities. 

Definition 1.1. Five functions (w £ , 9 e , p £ , q e , it 6 ) are called a generalized 
solution of model (NB) £ if they satisfy the regularity conditions in the 
domain Qt = x (0, T) 

w £ , D(s, w £ ), div x w £ , q £ , p £ , n £ , 6 £ , VJ £ G L 2 (Q T ), (1.4) 

at 

boundary conditions (10.91) in the trace sense, equations 

p £ + x £ a P div x w £ = 0, (1.6) 
tt £ + (1 - x £ )(a n dw x w £ - a es 6 £ ) = 0, (1.7) 

a.e. in Qt, and the integral equalities 

jf (a rP W ■ ^ - X^D^, ™ £ ) : D(x, ^) - p £ F ■ if + 
{(1 - x £ )a\B)(x, w £ ) - (q £ + ix £ )I} : D(x, (pj)dxdt + 
/ a T p £ (w £ ■ -£\t=o ~ v £ ■ <f\ t= ^jdx = (1.8) 
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for all smooth ip = ip(x,t) such that ip\on = <p\t=T = d<p/dt\ t=T = and 
jT ({a T c £ p 6 £ + ^div x i(; £ )^ - a%V x 9 E ■ + 
+ / (a r c£0§ + a £ e dw x w £ )£\ t=0 )dx = (1.9) 

for all smooth £ = £,(x,t) such that £|an = £|t=T = 0. 

In (11.81) by A : B we denote the convolution (or, equivalently, the inner 
tensor product) of two second-rank tensors along the both indexes, i.e., A : 
B = tr(B*oA) = J2l j=l A ij B ji . 

The first result on existence and uniqueness of solutions to model (NB) £ 
is the following lemma ([11]): 



Lemma 1.1. Whenever sequences {y/ct\(l — x £ )Vii>q}, {,/b~v }, {y/o~8^}, 
{y/a^(l — x £ )div x w £ }, {y/ct^x 6 div x wl} are uniformly bounded in L 2 (Q) with 
respect to e and F, ty, G L 2 (Q T ), model (NB) £ has a unique generalized 
solution (w € , 8 e , p e , q e , n e ) in the sense of Definition \l.lia nd 



m^(y/h~\\div x w £ {t)\\ 2Ml + ^fa{\\ W x w £ {t) || 2) n| 
dw £ 

+V^T\\-Qf( t )h,n + V^\\div x w £ (t)\\ 2) ne + Vc^\\6 £ (t)\\ 2 ,n) 
+V^f\\x £ V x 6 £ \\ 2 ,n T + s/a^ s \\(l - X £ )V x 6 £ \\ 2 ,n T 

+ v^llx^Va,-^— \\ 2 ,q t + y/a^\\x e div x —- Haor < (1-10) 
ot at 



where Co is a constant independent of e. 

§2. Formulation of the main results 

Suppose additionally that there exist limits (finite or infinite) 
lima p (e) = p*, lirna^e) = u , lunate) = x 0/ , lima^e) = x 0s , 

e\0 e\0 e\0 e\0 



OL f 

lim a v (e) = r] , lim a ds {e) = /3 Qs , lim a e /(e) = f3 0f , lim — = x xf 

e\0 e\0 e\0 e\0 OL^ 

hm — = /xi, Inn = p u hm = #u, Inn = p lf , 

e\0 E 2 e\0 Ct^ e\0 e\0 Cfy 



2 2 2 

hm = Ai, lim — ■ — = 771, hm = x ls , hm = x 2s , 

e\0 « M e\0 a M e\0 e\0 Ot\ 

V a V v a P v a 8f n v a 8s a 

hm— =7)2, hm— = p 2 , hm — - = piu hm — = p 2s . 

e\0 tt A e\0 a\ e\0 «a £ \° «A 

We assume that 

Assumption 2.1. Dimensionless parameters in the model (NB) £ satisfy 
restrictions 

0<T + /I 1 , X 0s , Xi /; A , P*, T]q] 

In what follows all parameters may take all permitted values. For exam- 
ple, if t = or p 0s = 0, then all terms in hnal equations containing these 
parameters disappear. 

The following Theorems 12.1142^1 are the main results of the paper. 

Theorem 2.1. Assume that conditions of lemma \l.l\ hold, dF/dt G L 2 (Qt), 
and that {w £ , 6 s , q £ , p £ , n £ } is a generalized solution in Model (NB) £ . 
The following assertions hold true: 

(I) U 

P*, Vo, A)/, Pos < 00, 

then 

max |||^(t)|,ya- X £ |V^ £ (t)|,(l-x £ )|V^ £ (t)||| 2 ,n<a, (2.1) 
||0 e |kn T + V^fWx £ V x 6 £ h,n T + ||(1 - X e )V,r || 2 ,n T < C , (2.2) 
\\q e \\2,n T + \\p e \\2,n T H H~aHkfi T + Ipltarv < Co- (2.3) 

(II) If 

A = x 0s = /ii = oo, < Ai, pi, ?7i, x ls < oo, 

and sequences {e^ 1 y /a x ~a^(l - x £ )^ w o}, fc -1 y/a^ v o} , \Jo^a^9 e Q } , 
{e' 1 yja^a^il — x £ )div x w £ ] }, {e^ 1 ^/a^a^x 6 div x w £ } are uniformly bounded 
with respect to e inL 2 (Q), then estimates (12.11) . (12.21) . (12. 3ft hold for pressures 
q £ ,p £ ,7T £ and re-normalized displacements and temperature 

w £ —> e~ 2 a^w £ , 6 £ — > e~ 2 a^8 £ , 
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with re-normalized parameters 



2 a ef 
<V 


otes - 


2 a 9s 


a v - 


2 a V 

^ 5 






2 ®-KS 


Otp 


, c 2 a P 




1 



2 2 ^-A 2 2 

C^/i 



a 0/ 



(HI) If 

\ = rj = oo, /ii < oo. 

T/ien /or displacements w e and temperature 9 £ hold true estimates ( 12.11) . 
(ED and if 

p*, (3 0f < oo, (2.4) 
then for the pressures q £ and p £ in the liquid component hold true estimates 

(O. 

// instead of restriction ( 12.41) hold true conditions 

f)<t> 3 V 

0<p 2 ,r ?2 <oo; ^ = 0; F = V$, — , | — | G L 2 (fi T ), (2.5) 
together with 

Assumption 2.2. Sequences {a\(l — x £ )^ w o}i {\/ qT ^ qT \ v o] > {v^j-c^oL 
{^/o^oa(1 — x e )div a .tUQ}, {y/a^aAX^div^tyo} are uniformly bounded with re- 
spect to £ in L 2 (fi). 

then are valid estimates 

max ||(1 - x e )\V x {a x w e )\ + |x £ ^(a A ^ £ )| + (2.6) 



|(1 - x £ )div x {a x w £ )\\\ 2 v + 1|(1 - X E )V x (a x 6 e )\\ 2 . nT < C , 

V "A 

which imply estimates (\2.3\i if Pof,Po s < oo. 

In (12.ip - (l2.6p Co is a constant independent of the small parameter e . 

Note, that for the last case {Ao = r] = oo, //i < oo} we can get same 
estimates (I2.1l) -( l2~3l) and (12.61) . if instead restrictions (12.51) and assumption 
12.21 we assume restrictions ( 12. 4ft together with 
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Assumption 2.3. 

l)F = F £ (x, t)(l - X e ), * = * £ (M(1 - X £ ), V2 < oo; 

2) sequences {a x (l-X e )Vw £ }, {^/a T axV £ }, { y Ja T a x 9 £ ] }, {^fa^(l-x £ )&w x w 
and {^a p a\x £ div x wl} are uniformly bounded with respect to e in L 2 (Q) and 
sequences {F £ }, {dF £ /dt}, {^ £ } and {<9\l/ £ /<9t} are uniformly bounded with 
respect to e in L 2 (Qt)- 

Theorem 2.2. Assume that the hypotheses in Theorem \2. 1\ hold, and 

Ho = 0; A , x 0s , p*, ?7o, A)/, /3 0s < oo. (2.7) 

T/ien functions w £ and 9 £ admit an extension u £ and d £ respectively from 
Q £ sT = Q £ s x (0, T) into such that the sequences {u £ } and converge 
strongly in L 2 (Qt) and weakly in L 2 ((0, T); W / 2 1 (fi)) to the functions u and d 
respectively. At the same time, sequences {w £ }, {6 £ }, {p £ }, {q £ }, and {vr e } 
converge weakly in L 2 (Qt) to w, 9, p, q, and tc, respectively. 

The following assertions for these limiting functions hold true: 

(I) If fi\ = oo then w = u, 9 = $ and the weak limits u, d, p, q, and tt 
satisfy in Qt the initial- boundary value problem 

d 2 u 

Top-^ + V(q + 7r)-pF= (2.8) 

div x {\ A s : B(x, u) + B s (div x u - + B{q}, 

Vo 

(r a e„ + &(1 - m ))% - *t f - = divJB ■ W) + % (2.9) 

r] dt p* dt r] dt 

—n + C S : D(x, u) + a s Jdiv x u - — ti) + a\q = 0, (2.10) 
Vo Vo 

—p+—ir + div x u-(l-m)—tf = 0, (2.11) 
P* Vo Vo 

q = P + W^+Pof&, (2-12) 



where 



p = mpf + (1 - m)p s , dp = mc p/ + (1 - m)c ps , m = J x{y)dy- 
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The symmetric strictly positively defined constant fourth-rank tensor Aq ; 
constant matrices Cq,Bq Bf, strictly positively defined constant matrix B^ 
and constants a\ are defined below by Eqs. (1 5 . 3 5 j) - (I5.37P and (15.401) . 

Differential equations (I2.8I) - (I2.12I) are endowed with initial conditions at 
t = and x £ Q 

Ou 

(T + p 0s )($-$ ) = 0, t (u-Uo)=t (—-v o ) = 0; (2.13) 

at 

and boundary conditions 

$(x,t) = 0, u(x,t) = 0, x e S, t>0. (2.14) 

(II) If the porous space is disconnected, then w = u and strong and 
weak limits u, p, q, tt together with a weak limit 6* of the sequence {x £ 9 £ } 
satisfy in Qj< equations (12. 8p . (12. 10)) — (12. lip , equation 

q = P + voP; 1 ^ + Po f f , (2.15) 

and heat equation 

r«c pf — + (r oCps + -)(!- m)- ------ = 

div x (B 9 ■ W) + ^. (2.16) 

Here 6* is defined below by Eqs. (15^21 - flB^g]) and A s , C S , B S Q B{, B® a s a{ 
are the same as in part (I) of the present theorem. 

The problem is endowed with initial and boundary conditions H2. 13)) - 
fl2HD. 

(Ill) If \i\ < oo then strong and weak limits u, d, w* , d* , p, q and it of 
the sequences {u 6 }, {$ £ }, {x e w 6 }, {x £ @ £ }, {p £ }, {q £ } ^nd {it 6 } satisfy the 
initial-boundary value problem infix, consisting of the balance of momentum 
equation 

r o (Pf + P'i 1 ~ m )^r) + V (<? + 7r ) " P F = ( 2 - 17 ) 

div x {X A s : D(x, u) + B s (div x u - — ■&) + B{q}, 

Vo 

where Aq ; Bq and B{ are the same as in ( 12.81) . continuity equation ( 12.101) . 
continuity equation 

— p H 7r + div x w^ — (1 — m)^-d — (m — l)div x u, (2.18) 

P* Vo Vo 
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state equation (I2.15p . heat equation (12.161) and Darcy's law in the form 



dw? du 



m- 



r f d 2 u 

J B 1 ( f i 1 ,t-r)-(-V x q + p f F-r p f -^)(x,r)dr (2.19) 



dt dt 

in the case of tq > and p\ > 0, Darcy's law in the form 
dw? du 



+ B 2 (p 1 )-(-V x q + p f F) (2.20) 



dt dt 

in the case of Tq = and, finally, Darcy's law in the form 



dt 



du 1 r 

5 3-^T + (ml-B 3 )- / (-V s q + p f F)(x,T)dT (2.21) 

dt T pf Jo 



in the case of p\ = 0. The problem is supplemented by boundary and initial 
conditions (12. 131) -( 12. 14ft for the displacement u and temperature $ of the 
rigid component and by the boundary condition 

w f (x, t) ■ n{x) = 0, (x,t)eS = dQ, t>0, (2.22) 

for the displacement of the liquid component. In Eqs. (I2.19P - (12.211) n(x) 
is the unit normal vector to S at a point x G S, and matrices Bi(pi,t), 

Theorem 2.3. Assume that the hypotheses in Theorem \2. 1\ hold, and that 

A = oo. (2.23) 

(I) If 

and one of conditions (12. 4p or (12. 5p holds true, then sequences {x £ u> £ }, 
{x £ @ £ }, {p £ }, and {q £ } converge weakly in L 2 {VLt) to w-f , 9* , p, and q respec- 
tively. The functions w £ and 8 s admit an extension u e and d e respectively 
from Q e s x (0,T) into Qt such that the sequences {u £ } and {$ 6 } converge 
strongly in L 2 (Qt) and weakly in L 2 ((0, T); W^i^l)) to zero and respec- 
tively and 

V if r o > and p\ > 0, then functions , 9* , d, p and q solve the 
problem (F x ) which consist of the state equation (I2.15p . continuity equation 

1 

— p + div x w f = 0, (2.24) 
P* 
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heat equation 

T ^ c vf% + <*(1 - ~ = ^v x {B e ■ W) + *, (2.25) 

Darcy 's law in the form 



Ot 



mv (x)+ [ B x (n u t-T)-{-V x q + Pf F){x,T)dT; (2.26) 
Jo 



r o = an d Pi > 0, then functions w-f , 9* , d, p and q solve the 
problem (F 2 ), which consist of equations (I2.15p . (12.241) . (I2.25j) and Darcy 7 s 
law in the form 

-jj r = B 2 ( f A 1 )-(-V s q + p f F) ] (2.27) 

finally, 

3) if To > and p\ = 0, then functions , 9> , *&, p and q solve the 
problem (F3), which consist of equations ( 12. 151) ( 12.241) . ( 12.251) and Darcy 's 
law in the form 

= mv (x) + ml -S3 • / {-V x q + p f F){x,r)dT. (2.28) 

dt T p f J 

In Eqs. ( 12. 26[) -( 12.281) matrices Bi(pi,t), .82(^1)? and B% are the same as 
in Theorem \2.2\ and 9* is defined below by Eqs. (15.421) - (15.481) . 

Problems F1-F3 are endowed with boundary condition (12.22[) for the dis- 
placement in the liquid component and initial condition (12 . 1 3[) . where /3q s = 0, 
and boundary condition (12. 14j) for the temperature in the solid component. 

(II) If 

H 0s , A)/, Pos, Pi < OO 

and conditions (12.51) together with assumption \2.2\ hold true, then the se- 
quence {a\u 6 } converges strongly in L 2 (Qt) and weakly in L 2 ((0, T); W 2 (^)) 
to function u, the sequence {tt £ } converges weakly in L 2 {VLt) to the function it 
and sequence {d £ } converges strongly in L 2 (VIt) and weakly in L 2 ((0, T); W 2 (p.)) 
to zero. The limiting functions u and ir satisfy the boundary value problem 
in the domain Q 



= div x {A s : B(x, u) + B s Q (div x u) + B[q - (g + tt) • 1} + pF, (2.29) 



15 



—7i + C* : B(x, u) + a s (div x u) + a\q = 0, (2.30) 

where the function q is referred to as given. It is defined from the corre- 
sponding of Problems Fi~F 3 (the choice of the problem depends on r and 
Hi), where $ = 0. The symmetric strictly positively defined constant fourth- 
rank tensor Ag, matrices Cq, Bq and B{ and constants a s Q and a\ are defined 
below by formulas (15.351) - (15.371) . in which we have r/o = r] 2 and A = 1. 

This problem is endowed with the homogeneous boundary conditions for 
the displacement u of the solid component. 

(HI) U 

fhf, fo s , X2s, Al, p 2 , Tl2<OG 

and assumption \2.2\ hold true, then after re-normalization 

w £ = a\w £ , -& e = a\-d £ 
we arrive at the problem, considered in Theorem \2.2\ where 

f = 0, fi x = —, A = 1, P* = P2, Vo = V2, 

ftos = fas, Pof = 02 f, X f = X 2 /, XOs = >?2s- 

(IV) Let 

/ii = oo, < Ai, pi, 771, flu, p lf , h u < 00 
and conditions (II) of Theorem \2.1\ hold. Then re-normalizing by 

w £ = a^E~ 2 w e , & = a^ 2 9 £ 
we arrive at the assumptions of Theorem \2.2\ where 

f = 0, #i = l, A = Ai, p*=pi, fj = r}i, 

ftos = ftls, A)/ = ftlf, i<0/ = *£ 0s = Ki s . 

(V) If 

H 1 = \ 1 = oo, 0<pi, 771, p ls , (3i f , h 1s <oo, 

then the corresponding problem for displacements {a^e~ 2 w £ } and tempera- 
ture {a^E~ 2 9 £ } has been considered in parts (I)-(II) of the present theorem. 
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Theorem 2.4. Assume that conditions of lemma \l.l\ hold and 
A , xo s , P*, Vo, Pof, Pos < oo, < fi . 

Then weak limits w, 6, p, q and n of sequences {w e }, {0 £ }, {p £ }, {q £ } and 
{7r e } satisfy in the following initial-boundary value problem: 

d^iv duo 

Top— + V(g + tt) = div x (A 2 : B(x, — ) + A 3 : B(x, w) + (2.31) 

B A div x w + 5{0 + / (A 4 (t - r) : B(x, io(cc, r)) + 
Jo 

B 5 (t - r)div x w(x, r) + B e 2 (t - t)6(x, r))dr) + pF, 

q = p + ^+(3 0f m6, (2.32) 

i r f 

—p + mdiv x w = - (C 2 (t - r) : B(x,w(x,r)) + (2.33) 
P* Jo 

a 2 (t — r)div x w(x, r) + a^t — t)9(x, r))dr. 



— tt + (1 - m)(div x w - —6) = - f (a d Jt - r)9(x, r) + (2.34) 
7/o Vo Jo 

C 3 (t — t) : D(x, w{x, r)) + a 3 (t — T)div x w(x, t))g?t, 



(r oCp + _(1 - -__-__ = ^(S ■ V.*) + *. (2.35) 

T/ie problem is supplemented by corresponding initial and boundary condi- 
tions (J2HSD -dHHD for 6 andw. 

In (l2~3TD - fl2T35|) A 2; A 3 and A 4 - fourth-rank tensors, B 4 , B 5 , B e , B{, B%, 
C 2 and C3 - matrices and a 2 , 03, a\ and a 2 - scalars. The exact expressions 
for these objects are given below by formulas (17. 38j) - (17.461) . The matrix B® 
is strictly positively defined. 

If the porous space is connected then A 2 is strictly positively defined sym- 
metric tensor. 

If the porous space is disconnected, which is a case of isolated pores, then 
A 2 = and the system (12.311) degenerates into nonlocal non-isotropic Lam'e's 
system with strictly positively defined and symmetric tensor A 3 . 
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§3. Preliminaries 



3.1. Two-scale convergence. Justification of Theorems 12. 1H2.41 relies 
on systematic use of the method of two-scale convergence, which had been 
proposed by G. Nguetseng [12] and has been applied recently to a wide range 
of homogenization problems (see, for example, the survey [S]). 

Definition 3.1. A sequence {p £ } C L 2 (Qt) is said to be two-scale convergent 
to a limit (p G L 2 (Qt x Y) if and only if for any 1-periodic in y function 
a = o~(x, t, y) the limiting relation 

lim / p e (x,t)a (x,t, x/e) dxdt = I I p(x } t,y)a(x,t } y)dydxdt (3.1) 
Jn T Jq t Jy 

holds. 

Existence and main properties of weakly convergent sequences are estab- 
lished by the following fundamental theorem p^2l [9] : 

Theorem 3.1. (Nguetseng's theorem) 

1. Any bounded in L 2 (Q) sequence contains a subsequence, two-scale 
convergent to some limit (p G L 2 (Qt x 50- 

2. Let sequences {p 6 } and {eV^} be uniformly bounded in L 2 (Qt)- Then 
there exist a 1-periodic in y function <p = p(x,t,y) and a subsequence {(f 6 } 
such that <f, Vyip G L 2 {VLt x Y), and ip £ and eV x (p e two-scale converge to <p 
and Vyip, respectively. 

3. Let sequences {<p £ } and {V^} be bounded in L 2 (Q). Then there exist 
functions (p G L 2 {VLt) andip G L 2 (OrxF) and a subsequence from {p 6 } such 
that ip is 1-periodic in y, V y ip G L 2 (Qt x Y), and ip £ and V x p> 6 two-scale 
converge to p and V x p(x,t) + V y ip(x,t,y) , respectively. 

Corollary 3.1. Let o G L 2 (Y) and o~ £ (x) := oixje). Assume that a se- 
quence {p 6 } C L 2 (Qt) two-scale converges to p G L 2 (Qt x Y). Then the 
sequence a e p £ two-scale converges to op. 

3.2. An extension lemma. The typical difficulty in homogenization 
problems while passing to a limit in Model (NB) e as e \ arises because 
of the fact that the bounds on the gradient of displacement V x w e may be 
distinct in liquid and rigid phases. The classical approach in overcoming 
this difficulty consists of constructing of extension to the whole Q of the 
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displacement field denned merely on The following lemma is valid due 
to the well-known results from [Tj [S] . We formulate it in appropriate for us 
form: 

Lemma 3.1. Suppose that assumptions of Sec. 1.2 on geometry of periodic 
structure hold, ip £ G W^iP^) and ip £ = on S e s = dQ £ s PI dQ in the trace 
sense. Then there exists a function a £ G W\ (fl) such that its restriction on 
the sub-domain fl £ s coincide with ip £ , i.e., 

(1 - x £ (x))(a £ (x) - i/f(x)) = 0, x G Q, (3.2) 

and, moreover, the estimate 

h £ \\ 2 ,n < CWrhm, HV^Ha.n < C|| V^lb,^ (3.3) 

hold true, where the constant C depends only on geometry Y and does not 
depend on e. 

3.3. Priedrichs— Poincare's inequality in periodic structure. The 

following lemma was proved by L. Tartar in [TH Appendix]. It specifies 
Friedrichs-Poincare's inequality for e-periodic structure. 

Lemma 3.2. Suppose that assumptions on the geometry of Q'y hold true. 

o 

Then for any function <p <EW\ (Q £ j) the inequality 

I \p\ 2 dx < Ce 2 j \V x p\ 2 dx (3.4) 

holds true with some constant C , independent of e. 

3.4. Some notation. Further we denote 

1) (®)y = fy ®dy, ($) y/ = J Yf My, ($)y s = J Ya $dy. 

2) If a and b are two vectors then the matrix a <g> b is defined by the 
formula 

(a ® b) • c = a(b ■ c) 

for any vector c. 

3) If B and C are two matrices, then B ® C is a forth-rank tensor such 
that its convolution with any matrix A is defined by the formula 

(B ® C) : A = B(C : A) 
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4) By l l i we denote the 3 x 3-matrix with just one non-vanishing entry, 
which is equal to one and stands in the i-th row and the j-th column. 

5) We also introduce 



./" =i(ny + P'*) = i(e i ®c i + c i 



where (ei, 62,63) are the standard Cartesian basis vectors. 

§4. Proof of Theorem EH] 

4.1. Let 

Ao, P*, Vo, Pof, Pos < 00. 

If restriction tq > holds, then ( 12. Ill — ( 12. 3p follow from the lemma 11.11 

Estimation of w £ and 9 £ in the case tq = is not simple, and we outline 
it in more detail. As usual, we obtain the basic estimates if we multiply 
equations for w £ by dw E /dt, equation for 9 e by 9 s , sum the result and then 
integrate by parts all obtained terms. The only two terms F ■ dw £ /dt and 
\I/ • 9 s heed additional consideration here. First of all, on the strength of 
Lemma 13.11 we construct an extension u e of the function w £ from Vt £ s into 
such that u £ = uf in D, £ s , u £ G W 2 1 (&) and 

Klkn < cyv^iiz.n < -^||(i - x e )V^>yxW% >Q . 



After that we estimate || io e || 2,0 with the help of Friedrichs-Poincare's 
inequality in periodic structure (lemma T3.2p for the difference (u £ — w £ ): 



\w £ h,n < \\u £ \\ 2 ,n + \\u £ - w £ \\ 2<n < \\u £ \\ 2 ,n + Ce\\x £ V x {u £ - w £ )\\ 2 . 

< \\u £ \\ 2 ,n + Ce||V a ,u e || 2 ,n + C{ea^)\\x £ ^V x w £ \\ 2)Q 

C -i 
< -=||(1 -X e )V^xW £ hv + C{sav 2 )\\x e V^xW £ hv- 



n 



The same method we apply for the temperature 9 £ : there is an extension 
$ £ of the function 6 £ from tt e s into tt £ f such that $ £ = 9 s in Q, £ , $ £ G W£{n) 
and 

U E h,Q < C\\V x $ £ \\ 2 , n < -^=||(1 ~ X e )v^V^1l2,n, 
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||r|| 2 ,n< -^\\(l-xnV^VJ £ \\ 2 ,n + C(eaJ)\\x £ V^VJ £ \\2.u> 



Next we pass the derivative with respect to time from dw £ /dt to p £ F and 
bound all obtained new terms in a usual way with the help of Holder and 
Grownwall's inequalities. Thus we have got estimates like estimates fll.lOp . 
but with right-hand side independent of To (or independent of e). 

4.2. The proof of this part of the theorem is obvious, because the re- 
normalization reduces this case to the case of p,\ — 1 and r = 0, which 
has been already considered above. 

4.3. Let A = oo, p t < oo and conditions H2.4[) hold true. It is obvious that 
estimates (12.11) . ( 12. 21) and estimates ( 12.31) for the pressures q £ and p £ in the 
liquid component are still valid. 

If conditions H 2 . 5 1) together with assumption 12.21 hold true, then the de- 
sired estimates (I2.6P follow from the basic equations for a\w £ and a\9 £ in 
the same way as in the case of estimates (12.11) - (12. 3p . The main difference 
here is in the term p £ F ■ a\dw £ / dt, which now transforms to 



dw £ , w ^ „ dw 



£ 



T = p f F- a x —- + (p f - Pf )(l - X £ )F ■ a x 



dt ' yr] * ' A dt ' 
The integral of first term in T transforms as 

f* f „^ dw £ , , f* r dw £ , , 

Pf / \7<P ■ ax—^ — dxar = —pf / / <Paxdiv x — — axdr 



/"* f dw £ dw £ 

= -Pf I (X £ ■ $axdiv x — h (1 - x £ ) ■ ^a x div x ^—)dxdT 

Jo Jn or dr 

= -Pf I (X £ ■ ®axdiv x w £ + (1 - x £ ) • <5>a x div x u £ )dx+ 
Jn 

Pf (X £ ■ $> T axdiv x w £ + (1 - x £ ) ■ ^ T axdw x u £ )dxdr 

Jo Jn 

and is bounded with the help of terms 

/ (^Ka^XaAdiv^) 2 + (1 - x £ )\axV x u £ \ 2 )dx, 
Jn 

which appear in the basic identity after using the continuity equations. 
The integral of the second term in T is bounded with the help of the term 



"a 



/ (l-X E )|V B ti 6 | 2 da; 
Jn 
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in the same way as before. 

Estimates (12. 3p follow now from f)2.6p if /3 /,/?os < oo. 



□ 



§5. Proof of Theorem [Ml 

5.1. Weak and two-scale limits of sequences of displacement, 
temperatures and pressures. On the strength of Theorem 12.11 the se- 
quences {0 e }, {p e }, {q £ }, {tt 6 } and {w £ } are uniformly in e bounded in 
L 2 (Q T ). Hence there exist a subsequence of small parameters {e > 0} and 
functions 9, p, q, 7r and w such that 

6 s —> 6, p e — > p, q £ — »• q, ir £ — > n, w £ — w 

weakly in L 2 (Q T ) as e \ 0. 

Due to Lemma [3.11 there is a function u £ 6 L°°((0,T); W / 2 1 (r2)) such that 
w e = in fi s x (0, T), and the family {it e } is uniformly in £ bounded 
in L°°((0, T); W^fi)). Therefore it is possible to extract a subsequence of 
{£ > 0} such that 

w £ -> w weakly in L 2 ((0, T); W^O)) 

as £ \ 0. 

Applying again the same lemma 13.11 we conclude that there is a function 
$ £ G L 2 ((0,T); W^O)) such that ^ = fl £ in fi s x (0,T), and the family {$ £ } 
is uniformly in e bounded in L 2 ((0, T); iy 2 1 (f2)). Therefore it is possible to 
extract a subsequence of {e > 0} such that 

d £ -> t? weakly in L 2 ((0, T); W^ft)) 

as £ \ 0. 
Moreover, 

X £ « M D(^, ™ £ ) -> 0, x^/Vfl 6 -> (5.1) 

as £ \ 0. 

Relabelling if necessary, we assume that the sequences converge them- 
selves. 

On the strength of Nguetseng's theorem, there exist 1-periodic in y func- 
tions Q(x,t,y), P(x,t,y), U(x,t,y), Q(x,t,y), W(x,t,y), Q s (x,t,y) and 
U(x,t,y) such that the sequences {0 £ }, {p £ }, {tt £ }, {q 8 }, { w£ }, {V^ 6 } and 
{V x u £ } two-scale converge to Q(x,t,y), P(x,t,y), H(x,t,y), Q(x,t,y), 
W(x, t, y), + V y s (a;, t, y) and V x u + V y U(x, t, y), respectively. 
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Note that the sequence {div x w 6 } weakly converges to div x w and $,\u\ G 

o 

L 2 ((0, T); W\ (fi)). Last assertion for disconnected porous space follows from 

o 

inclusion $ £ , \u £ \ G L 2 ((0,T); W\ (fi)) and for the connected porous space it 
follows from the Friedrichs-Poincare's inequality for u £ and d £ in the e-layer 
of the boundary S and from convergence of sequences {u £ } and {i? 6 } to u 
and $ respectively strongly in L 2 (Qt) and weakly in L 2 ((0, T); W^fi)). 

5.2. Micro- and macroscopic equations I. 

Lemma 5.1. For all x £ Q and y G Y weak and two-scale limits of the 

sequences {0 £ }, {p 6 }, {tt £ } ; {q £ }, {w £ }, {V x $ £ } and {V x u £ } satisfy the 
relations 

Q = ixq, Q = P + ^dP/dt + p ofX Q; (5.2) 

ill + (1 - X )((div x u + <%17) - = 0; (5.3) 

dityW" = 0; (5.4) 

W = X W+(l-x)u; (5.5) 

e = x e + (i- x >?; (5.6) 

q = V + VtiP?dp/dt + P of f , 6f = (e) Yf ; (5.7) 
±n + (1 _ + (^C7) Ks - (1 - m)^tf = 0; (5.8) 

—■K + —p+dzv x w-(l-Tn)^^ = 0. (5.9) 

Proof. The weak and two-scale limiting passage in Eq. (II. 6p yield that Eq. (15.71) 
and the second equation in (15. 2p . 

In order to prove first equation in (15. 2p . into Eq. fll.8p insert a test function 
i/> £ = (x, t,x/e), where ij)(x, t, y) is an arbitrary 1-periodic and finite on 
Yf function in y. Passing to the limit as e \ 0, we get 

V y Q(x,t,y) = 0, yeY f . (5.10) 

Next, fulfilling the two-scale limiting passage in the equality 

(1 - X £ )P £ = 

we arrive at 

(1- X )P = 
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which along with Eqs. fl5.10p and second equation in (J572]) justifies first equa- 
tion in Q . 

Eqs. ( 15.31) . (15. 4p . (15. 8p . and (15.91) appear as the results of two-scale limiting 
passages in Eqs. fll.6p and (II. 7p with the proper test functions being involved. 
Thus, for example, Eq. (15.81) is just a subsequence of Eq. ()5.3p and Eq. (l5.9p is 
a result of two-scale convergence in the sum of Eq. (ll.6p and Eq. (ll.7p with 
the test functions independent of the "fast" variable x/e. Eq. (15.41) is derived 
quite similarly if we represent Eq. (11.61) in the form 

— p £ + div x w £ = (1 - X £ )div x u £ , (5.11) 

Otp 

multiply by an arbitrary function x\) £ = exp (x,t,x/e) and then pass to the 
limit as e \ 0. 

In order to prove Eqs. (l5.5l )- (l5.6p it is sufficient to consider the two-scale 
limiting relations in 

(1 - x £ )(™ £ - u £ ) = 0, (1 - X £ W - #) = 0. 

□ 

Lemma 5.2. For all (x,t) G Qt the relations 

\ A y U = V y U, yeY s , (5.12) 

(\ Q B(y, U) - n • I + A D(x, u) + —q • l) ■ n = 0, j/67, (5.13) 

m 

hold true. Here n is a unit normal to 7. 

Proof. Substituting a test function of the form if) £ = e\\) (x,t,x/s), where 
ip{x, t, y) is an arbitrary 1-periodic in y function vanishing on the boundary 
dfl, into Eq. (11.81) . and passing to the limit as e \ 0, we arrive at the following 
microscopic relation on the cell Y: 

div y {A (l - x)(D(y, U) + D(x, «)) - (n + -q X ) • 1} = 0, (5.14) 

m 

which is clearly equivalent to Eqs. (15. 121) and (1 5 . 1 3 [) in view of Eqs. (l5.2p . □ 

In the same way using additionally continuity equations ( 11.61) and ( 11.71) 
one get from Eqs. ( 11.91) 
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Lemma 5.3. For all (x, t) G Qt the relations 

A y Q s = 0, yeY s , (5.15) 
— = -y^ • n, ye 7 (5.16) 

Now we pass to the macroscopic equations for the solid displacements. 

Lemma 5.4. Let p = mpf + (1 — m)p s , = (W)y } - Then functions 

u,w* , q,7r,9 : ' ,d satisfies in Qt the system of macroscopic equations 

T oPf-Q^ + T oPs(~L-m)—-pF= (5.17) 
div x {X ((l - m)B(x, u) + (3{y, U)) Ys ) - (q + tt) • I}, 

roe pf ^ + (r c ps + -)(!- m)- ------= (5.18) 

x 0s div x {(l - m)V x d + (V y 6 s )yJ + 

Proof. Eqs. fl5T7D and fl5TT8D arise as the limit of Eqs.([LHD and CCHD with 
test functions being finite in VLt and independent of e. In Eq. fll.9l) we have 
used continuity equations (\1.6\i and (\1.7\i . □ 

5.3. Micro- and macroscopic equations II. 

Lemma 5.5. // p\ = oo then u = w and 6 = d. 

Proof. In order to verify, it is sufficient to consider the differences (u £ — w £ ) 
and (6 s — $ e ) and apply Friedrichs-Poincar'e's inequality, just like in the 
proof of Theorem 12.11 □ 

Lemma 5.6. Let pi < oo and V = xdW/dt. Then 

dV 

T 0Pf —-p f F = p 1 A y V-\7 y R-V x q, yeY f , (5.19) 

dO 3 dP 

^oc p/ ^ = x lyUl A s/ e + ^— + vl/, y G Yf, (5.20) 

V = —, e = #, ye 1 (5.21) 
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in the case [i\ > 0, and 
dV 

T 0Pf — = -V y R-V x q + p f F, yeY f , (5.22) 

(xVT - w) - n = 0, y G 7 (5.24) 

in the case \i\ = 0. 

In .Eg. ( 15.24)) n zs the unit normal to 7. 

Proof. Differential equations (15.191) and (15.221) follow as e \ from inte- 
gral equality ( 11. 8ft with the test function if) = ip(xe~ l ) ■ h(x,t), where (p is 
solenoidal and finite in Yf. 

The same arguments apply for the Eq. ( 15.201) and Eq. (l5.23p . The only one 
difference here is that we use the continuity equation (II. 6p to exclude the 
term x £ d.iv x (dw e /dt). 

First boundary condition in (15.211) is the consequence of the two-scale 
i_ 1 

convergence of {a^V x w £ } to the function \i\ V y W(x, t, y). On the strength 
of this convergence, the function V y W(x, t, y) is L 2 -integrable in Y. As 
above we apply the same argument to the second boundary condition in 
( I5.2ip . The boundary conditions (I5.24p follow from Eq. (15.4)) . □ 

Lemma 5.7. If the porous space is disconnected, which is the case of isolated 
pores, then u = w. 

Proof. Indeed, in the case < jii < 00 the systems of equations (15. 4p . 
( 15TT9|) . and (15T2T|) . or Q, fl222D, and (1041 have the unique solution V = 
du/dt. □ 

5.4. Homogenized equations I. 

Lemma 5.8. If fii = 00 then w = u, 9 = d and the weak limits u, d, p, q, 
and tt satisfy in the initial-boundary value problem 

ToP-Q^ + V(q + 7c)-pF= (5.25) 
div x {\ A^ : D(x, u) + B s Q (div x ii - + B{q), 
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(r„c t + &(1 - m))% - Sm.% - *4 = • V*) + •, (5.26) 

— + C S : D(x, u) + agCdi^t* - — ■&) + a?g = 0, (5.27) 
Vo ?7o 

— p+— ir + div x u - (1 - m)— -9 = 0, (5.28) 
P* ?7o Vo 

dp 

q = p + v oP* ~^ + Pofmfi, (5-29) 

where the symmetric strictly positively defined constant fourth-rank tensor 
Aq, constant matrices Cq, Bq B{, strictly positively defined constant matrix 
B® and constants a s Q a\ are defined below by formulas (1 5 . 3 5 [) - (I5.37p . (15.401) . 

Differential equations (I5.25l) - (l5.26p are endowed with initial conditions at 
t = and x e Vt 

Ou 

(r + faW - #„) = 0, r (u - u ) = r (— - v ) = 0; (5.30) 
and boundary conditions 

0(x,t) = O, u(x,t) = 0, xeS, t>0. (5.31) 
Proof. In the first place let us notice that u = w and 9 = •& due to Lemma 



The homogenized equations (15.251) follow from the macroscopic equations 
(I5.17p . after we insert in them the expression 

(B(y, U)) Ys = A? : D(x, u) + B s (div x u - + B{q. 

Vo 

In turn, this expression follows by virtue of solutions of Eqs.f l5.3p and (15.141) 
on the pattern cell Y s . Indeed, setting 

3 

U=J2 U ij (y)D ij + U (y)(div x u - ^0) + -U 1 (y)q 
tftx Vo 

and 

3 

n = A J2 U i ^(y)D ij + Tl (y)(div x u - ^0) + -U 1 (y)q 1 
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where 

ldu 1+ du J 



2 ' &Xj ' i).V, 



we arrive at the following periodic-boundary value problems in Y: 

div„{(i - x)(P(y, u ij ) + J lJ ) - n ij ■ i} = o, 

^IF + (1 - x)div,J7« = 0; 



(5.32) 



(5.33) 
(5.34) 



dw y {\o(i - xWv, u ) - n • 1} = 0, 

in o + (l-x)(div y C/o + l)=0; 

di Vj/ {A (i - xMv, Ui) - (n x + x) ■ 1} = o, 

^ + (1 - X )div 2/ t/ 1 ) = 0. 

On the strength of the assumptions on the geometry of the pattern "liq- 
uid" cell Y s , problems (15.32)) — ( 15.341) have unique solution, up to an arbitrary 
constant vector. In order to discard the arbitrary constant vectors we de- 
mand 

(U ij )y s = (U ) Ys = (U 1 ) Ya = 0. 

Thus 



3 3 

\ s — 11 

^0 



J2j ij ®J ij + K H = J2((l-x)D(y,U i i)) Y ®J ij . (5.35) 

i,j=l i,j=l 



Symmetry and strict positiveness of the tensor Aq have been proved in [TO] . 

Finally, Eqs. (l5"^7D - fl5^9D for the pressures follow from Eqs.j52j)- CT) 
and 

B s = (B(y,U )) Ya , B^-iBiy,^))^, a{ = -(diVyUx)^, (5.36) 

3 

a% = 1 - m + (dw y U ) Ys , C s = (diy y U ij ) Ys J». (5.37) 

Now for i = 1,2,3 we consider the model problems 

A„e? = 0, ye Y s , (5.38) 

dot 

— — = -e r n, yG7 
on 
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and put 



3 



i=i 



Then G s solves the problem (15.150 — (15.160 and 



3 



B e = x 0s ((l - m)l + ^(Vj,8?) n <g> e;). (5.40) 

All properties of the matrix 5 are well known ( see [H], [8]). □ 

Lemma 5.9. // the porous space is disconnected, then w = u and the weak 
limits Q* , u, p, q, and tc satisfy in fl T equations (15.251) . (15.271) - (15.281) . 
(15.71) . where Aq ; Cq,Bq B{, B® , a s Q and a{ are the same as in Lemma \57R 
and heat equation 

r«c pf — + (r c ps + -)(!- m)- ------ = 

div x (B e ■ W) + ^. (5.41) 
If Pof = 0; then for the case fii > 0, r > 

6 f (x,t) = md{x,t) + f b e f {t-T)h(x,T)dT, (5.42) 

Jo 

where 

h(x,t) = V(x,t) - T c pf —(x,t); 

for for the case r = 

6 f (x,t)=m$(x,t)-c e f y(x,t); (5.43) 
and, finally, for the case /ii = 

171 

9 f (x,t) =m& (x) + / V(x,T)dT. (5.44) 

T 0Cpf Jo 

Here b 9 j-{t) and c® are defined below by formulas (I5.45p . (15.461) . 

If v o — 0, then Q* is defined by the same formulas (15.421) -( I5.46P where 
T c p f is replaced by (r Q c p f + -^) and \1/ is replaced by (\1/ + ^q). 

If > and /3qj > 0, then 9* is defined by (15.421) where h(x,t) and b®(t) 
are defined below by formulas ( 15.471) and ( 15.481) . 



The problem is endowed with initial and boundary conditions ( I5.30fl 
(J53U). 
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Proof. The only one difference here with the previous lemma is the heat equa- 
tion for because 8 ^ The function 9? now is defined from microscopic 
equations O, (TOUD and (IB^TD . 

If /3 0/ = 0, then the system ([53]) . fl5T20l) and f[5T2Tj) is decoupled and 0' 
is represented by (I5.42I) - (I5.44I) where 

&}(*) = <©{>!>, 4=<©oV (5.45) 
and 6{, 6q are solutions of periodic in y problems 

qqS 

T o c pf-gf = xi/iiA^ef, y e Y f , (5.46) 

e{(i/,0) = l, yey /; ef = o, ye 7 , 

and 

K^nAyBl = 1, 2/ G 17; 6^ = 0, ?/ G 7. 

The pressure p is defined from (15. 7p . 

If > and /3 / > 0, then we look for the solution of (15.21) . (15.201) and 
(I5.2ip in the form 

G(x,t,y) =ti(x,t)+ [ e f (y,t-T)h(x,r)dT, 

Jo 

P(x,t,y) = h {x,t)+ [ P f {y,t-T)h{x,T)dT. 

Jo 

Functions h(x,t) and ho{x,t) solve a system 



Pof,dh , d'd 

— (-QT + h )-T0C p fh = T C pf --% 

p* ot 111 
and 1-periodic in y functions Qf(y,t) and Pf(y,t) solve a problem 

r c pf — = > (m A y ef + — — , yEY f , 
p* dt 
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(5.47) 



r (e'-l)| t= o = (P-l)|t = o = 0, yeY f ; /Xi6> = 0, ye 1 . 
This problem has a unique solution due to estimates 

f [ { Hlfil \vef(y,T)\ 2 + ^\^(y,r)\ 2 }dydr+ 

JO JYj P* OZ 

r 1 1 m 1 

/ {-T c pf \ef(y,t)f + —\pf(y,mdy = -(r c p/ + -). 

JY f ^ L P* L P* 

Thus, Qf is defined by (15 .421) where h(x,t) is defined by (15.471) and 

6}(f) = (&%. (5.48) 

□ 

5.5. Homogenized equations II. 

Let fi\ < oo. In the same manner as above, we verify that the weak limit 
u of the sequence {u 6 } satisfies some initial-boundary value problem likes 
problem (15.251) - (15.31!) because, in general, the weak limit w of the sequence 
{w £ } differs from u. More precisely, the following statement is true. 

Lemma 5.10. If fii < oo then the weak limits u, , 6? , p, q, and n of 

the sequences {u 6 }, {x £ w £ }, {x £ 9 £ }, {$ £ }, {p £ }, {q £ }, an d {^ £ } satisfy the 
initial-boundary value problem in Vtx, consisting of the balance of momentum 
equation 

^Pf-^T + P'O- ~ m )^) + V (<? + *)-pF= (5.49) 
div x {X Al : B(x, u) + B s Q div x u + B{q), 

where Aq, Bq and B{ are the same as in (15.251) . the continuity equations 
(I5.27p . equation 

— p H 7r + div x w-f + -^-'d = (m — l)div x u, (5.50) 

P* Vo Vo 

equation (15.71) . heat equation ( I5.4ip and Darcy's law in the form 
i)w f du r ' > y2 " 



dt dt 



r d 2 u 

J B 1 ( f i 1 ,t-r)-(-V x q + p f F - To p f —)(x, T )dT (5.51) 
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in the case of Tq > and fix > 0, Darcy's law in the form 



dw? du 



+ J B 2 (/i 1 ) • (-V x q + Pf F) 



(5.52) 



dt dt 



in the case of r = and, finally, Darcy's law in the form 




1 



(mI-B 3 )- [ (-V x q + p f F)(x,r)dr (5.53) 



dt 



in the case of pi = 0. The problem is supplemented by boundary and initial 
conditions (15.301) ( T5.31H for the displacement u and temperature d of the 
rigid component and by the boundary condition 



for the displacement w-f of the liquid component. In Eqs. (15. 51 p - (15.541) n(x) 
is the unit normal vector to S at a point x e S, and matrices Bi(fii,t), 
-^2(^1); and B 3 are given in fWf . 

The proof of these statements repeats the proof of Lemma 5.8 in [TU] . 



6.1. Weak and two-scale limits of sequences of displacements, 
temperatures and pressures. 



X 0s , (3 f, A)s, yUl < OO 

and one of the conditions (12 ,4p or (I2.5P holds true. Then on the strength 
of Theorems 12.11 and 13.11 we conclude that sequences {x £ w £ }, {x £ 9 £ }, {p £ } 
and {q 6 } two-scale converge to x{y)W(x, t, y), x(y)®(x, t, y), P(x, t, y) and 
Q(x,t,y) and weakly converge in L 2 {VLt) to w-f , 9% p and q respectively. A 
sequence {u 6 }, where u £ (x, t) is an extension of w £ (x, t) from the domain Q £ 
into domain Q, strongly converges in L 2 (Qt) and weakly in L 2 ((0, T); (^)) 
to zero and a sequence {^ e }, where i9 £ (x,t) is an extension of 8 £ (x,t) from 
the domain Q £ s into domain Q, strongly converges in L 2 (Qt) and weakly in 
L 2 ((0, T); iy 2 1 (f2)) to Moreover a sequence {V x $ £ (x, t)} two-scale converge 
to X7 x $(x,t) + V y e s (x,t,y). 



w 



f (x, t) ■ n(x) = 0, (x,t)eS = dn, t>0, 



(5.54) 



§6. Proof of Theorem EH 



I. Let 



II. If 



XQs, A)/, fas, yUl < OO 
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and conditions (12. 5p together with assumption 12.21 hold true, then due to 
estimates (I2.3P and (12.61) the sequences {ot\U £ } and {ft 6 } converge strongly in 
L 2 {VL T ) and weakly in L 2 ((0, T); W^^fl)) to a function u and zero respectively, 
and the sequence {tt £ } converges weakly in L 2 (Qt) to a function 7r. 

III. If 

x 2s , x^ 1 , (3 2f , p 2 s, Ai, p 2 , r] 2 <oo 

and assumption 12.21 holds true, then due to estimates (I2.1I) - (I2.3I) for the re- 
normalized displacements and temperature sequences {x e C(\w £ }, {x e a\0 £ }, 
{p £ }i {q £ } an d { n£ } weakly converge in L 2 (VLt) to , Qf , p, g and ff respec- 
tively and sequences {aAit e } and {cka 1 ^} converge strongly in L 2 (Qt) and 
weakly in L 2 ((0,T); W / 2 1 (f2)) to functions u and $ respectively. 

IV. If /i! = cxd and 

< Ai, pt, 77i, /3i s , /?!/, Xi a < 00, 

then on the strength of part 2) of Theorem 12. II the sequences {a^e~ 2 x £ i^ 6 }, 
{a^e^x 6 ^ 6 }, {p 6 }, {q 6 } and {7r e } weakly converge in L 2 (VL T ) to , 6^, p, 
g and 7T respectively, and the sequences {a^e~ 2 u 6 } and {a^ 2 d £ } converge 
strongly in L 2 (fl T ) and weakly in L 2 ((0, T); W / 2 1 (r2)) to functions u and $ 
respectively. 

As before in previous section §5, we conclude that 

o 

tf, |t*|,0, |u| e L 2 ((0,T); Wj 1 (fi)). 
6. 2. Homogenized equations. 

I. We construct a closed system of equations for the velocity dw^ jdt in 
the liquid component pressures p and q as before ( proof of Theoren i2.2p . 
consisting of state equation (12.151) continuity equation (I2.24p . heat equa- 
tion f)2.25p . modifications of Darcy's law (I5.5ip - (I5.53P ( in which we have 
u(x,t) = 0), boundary condition (15.541) for dw-f/dt and boundary and ini- 
tial conditions (I5.30p . (I5.3ip for the temperature of the solid component. 
To derive continuity equation (I2.24p we use equation (15. lip and to derive a 
homogenized heat equation (I2.25P we rewrite the heat equation (11.91) using 
continuity equation (11.61) . Note that in the previous section to get the same 
homogenized equation for the temperature of the solid component we have 
used both continuity equations (I1.6P and (II. 7p . Here we do not express the 
term (1 — x £ )div(du e /dt), because it converges to zero. 

We entitle the above described systems as Problem Fx, F 2 , or F 3 depend- 
ing on the relations between parameters r and Hi. 
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1) t > 0, n\ > for the problem Fi, 

2) t = for the problem F 2 ; 

3) Hi = for the problem F 3 . 

II. In this situation we observe that the limiting displacements in the rigid 
skeleton are equal to zero. In order to find a more accurate asymptotic of the 
solution of the original model, we use again the re-normalization. Namely, 
let 

ur — > a\W e . 

Then new displacements satisfy the same problem as displacements before 
re-normalization, but with new parameters 

a v — > a^a^ 1 , ot\ — > 1, a T — > a T a^. 

Thus we arrive at the assumptions of Theorem 12.21 the limiting functions 
u(x,t), 7i(x,t), H(x,t,y), and U(x,t,y) satisfy the system of micro- and 
macroscopic equations (15. 3p . (I5.14p . (15.81) and (15.171) . in which the pressure 
q is given by virtue of one of Problems Fy- F 3 , and /3 0s = 0, r = 0, r] = r] 2 , 
A = l. 

Hence for u(x, t) and ir(x, t) there hold true the homogenized momentum 
equation ( 12.291) . homogenized continuity equation ( 12.301) and a homogeneous 
boundary condition for u. □ 

§7. Proof of Theorem EH 

7.1. Weak and two-scale limits of sequences of displacement, 
temperature and pressures. 

On the strength of Theorem 12.11 the sequences {w £ }, 9 £ , {p £ }, {q e } and 
{7r £ } are uniformly in e bounded in L 2 (Qt)- Then there exist a subsequence 
from {e > 0} and functions p, n, q, and w such that as e \ 

w £ ^w, e e -> 9, f ^p, q £ -*q, ir £ -> tt L 2 (VL t ). (7.1) 

Moreover, since A ,/io > then the bounds (12. lj) — (12. 3p imply 

V x w £ — ► X7 x w, V x 9 £ — ► V x 9 weakly in L 2 (Vt T ). (7.2) 

e\0 e\0 

Due to limiting relations (17. ip . (17. 2p and Ngutseng's theorem, there exist 
one more subsequence from {e > 0} and 1-periodic in y functions P(x, t, y), 
H(x, t, y), Q(x,t,y), Q(x,t,y) and U(x,t,y) such that the sequences {p £ }, 
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{7r £ }, {q £ }, {V9 £ } and {Vw £ } two-scale converge as e \ respectively to 
P, IT, Q, V x 9 + V^e , and V x w + VyU. 

7.2. Micro- and macroscopic equations. 

Lemma 7.1. Two-scale limits of the sequences {p £ }, {^ £ }, {q £ }, {V8 £ } and 
{Vw £ } satisfy inYr = Y x (0, T) the following relations 

— n + (1 - x)(div x w + div v W - —&) = 0; (7.3) 
Vo Vo 

—P + x{div x w + diVyW) = 0; (7.4) 

P* 

dwy(xiio(B(x, + D(y, — )) + (1 - x)A (O(x, to) + D(y, W))) (7.6) 

-V„(Q + n) = 0; 

divy( X x of (v x e + v„e) + (i - x)x 0s (v^ + v„e)) = o. (7.7) 

Lemma 7.2. TTie weak limits p, n, q w satisfy in Qt the following system 
of macroscopic equations: 

— 7T + (1 - m)div x w + (div y W) Ys - (1 - m)—6 = 0; (7.8) 
Vo Vo 

— p + mdiv x w + (dm. WW. = 0; (7.9) 

P* 

q = p+-% + (3 of m0; (7.10) 
p* ot 

O'^xv duo 
T *P-QV + V( - q + 7r )-P F = dw x (no(rriB(x, —) + (7.11) 

3W 

-Qf))y } ) + V(l " m)D(x, w) + (D(y, W))y,)); 

w+ff 1 -"))!-^!-!!-- <" 2 » 

div x (x f(mV x 6+ (V y G) Yf ) + *b*((l - m)V x + (V w e)yJ). 
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Proofs of these statements are the same as in lemmas 15.11 - 15.41 
7.3. Homogenized equations. 

Lemma 7.3. Weak limits p, n, q, 9 and w satisfy in Qt the following system 
of homogenized equations: 

Top— + V(q + tt) = div x (A 2 : B(x, — ) + A 3 : B(x, w) + (7.13) 

B 4 div x w + B°6+ [ (A 4 (t-r) : B(x,w(x,t)) + 
Jo 

B B (t - r)div x w(x, t) + B 2 (t - r)9(x, r))dr) + pF, 

q = P+~ + Pofm9, (7.14) 

1 /"* 

—p + mdiv x w = - (C 2 (t - r) : B(x,w(x,t)) + (7.15) 
P* Jo 

&2(t — r)div x w(x, t) + a x (t — t)9(x, r))dr. 

— 7i + (1 - m)(div x w - ^9) = - [ (a e 2 {t - t)6{x, t) + (7.16) 
7/o Vo Jo 

C 3 (t — t) : B)(x, w(x, r)) + a 3 (t — r)div x w(x, r))dr, 



(r oCp + -(1 - m))- --!_-__ = ^(S ■ + *. (7.17) 

i/ere A2, A3, A4 - fourth-rank tensors, B4, B^, Bq, B\, B 2 , Ci, C3 - 
matrices and a 2 , a 3 , a\, a 2 - scalars. The exact expressions for these objects 
are given below by formulas ( I7.38p - (17.461) . 



Proof. Let 



dw 

Z(x, t) = fj, Q B{x, —) - A D(a;, w), Z {j = e t ■ (Z ■ ej) 



z(x, t) = ^ Zi(x, t)ei = (x 0/ - x f)V x 9, z {x, t) = div x w. 

i=l 
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As usual we look for the solution of the system of microscopic equations 
dZ3D-(I72D in the form 

W [ t [W (y,t-T)z (x,r) + W 9 (y,t-r)e(x,r) (7.18) 
Jo 

3 

+ Y / W ij (y,t-r)Z ij (x,r)]dT, 

P = X [ [P°(y,t-r)z (x,r)+ (7.19) 
Jo 

3 

+P e (y,t - t)0(x,t) + P i3 (v^- r)Z i3 {x,r)\dr, 

3 

Q = x(Qo(y) ■ zo(x,t) + Q j i{y)Z ij {x 1 t) (7.20) 

+Q e (y)6(x,t) + [ [Q (y,t-T)z (x,r)+Q e (y,t-r)e(x,r) 
Jo 

3 

+ J^ / Q ii (y,t-T)Z ij (x,T)]dT), 

n=(l- X ) [ t [J2n ij (y,t-T)Z ij (x,T)+ (7.21) 
U°(y, t - t)z (x, t) + n e (y, t - t)9(x, t)] dr, 

3 

9 = ^6%)^,*), (7.22) 



i=l 



where 1-periodic in y functions W°, W\ W ij , P°, P e , P", Q , Qg, g?', Q°, 
<5 e , Q u , n°, IT 61 , IP- 7 , l satisfy the following periodic initial-boundary value 
problems in the elementary cell Y: 
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Problem (7) 

dW ij 

(1 - X )(A D(y, W) - ((1 - X )IF + XQ 13 )1)) = 0; (7.23) 

-P ij + xdiv^W^' = 0, Q ij = P ij + 

p* at 

In* + (i _ x )div y W^ = 0, W*(f/, 0) = W?'(i/); (7.24) 

di Vy ( X (/ioO(y, Wf) + J iJ - gj'l)) = 0, (7.25) 

x(Q« + ^div,^') = 0. (7.26) 
Problem (77) 

div,(x(^ ©(y,^ r ) + 

(1 - X )(AoD(y, W°) - ((1 - X )n° + xQ°)I)) = 0; (7.27) 

x(-P° + dW y W + l) = 0, Q° = P° + ^^; (7.28) 

(1 - x) (— n° + div y W° + 1) = 0; (7.29) 

Ty°(y,0) = W°(y), div,( X (//oD(y,W-°)-QoI))=0, (7.30) 

X (Qo + ^o(div 2/ W° + l)) = 0. (7.31) 

Problem (777) 

div,(x(/i D(y,^ r ) + 

(1 - X )(AoD(y, W e ) - ((1 - X )n e + X Q e )I)) = 0; (7.32) 

x (-P» + dxv y W) = 0, Q° = P e + + PofX; (7.33) 

J)* Ob 

(1 - x) (-n e + divj, W e - ^) = 0; (7.34) 

W e (y,0) = W e (y), div y (x(iM&(y,W o )-Q e o I))=O, (7.35) 

X (Q e + v div y W e -[3 0f )=0. (7.36) 

Problem (IV) 

div y ( X H 0f + (1 - X )xo s )V,0 J + X e t ) = 0. (7.37) 
38 



Then 
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A 2 = Horn J ij g> J ij + n A f Q , 

3 

A£ =HoJ2 (°(^ ® J ^ ( 7 - 38 ) 



A (l - m) J] J ij <g> J ij - A A£ + /x A{(0), 



ij'=l 

A 4 (t) = Voj t M(t) - A A{(t); (7.39) 

A{(t) = ^{^Wy.-^(w,0)>y/ + Ao<D(y,W«(y,t))>y.}® J«; (7.40) 

B 5 (t) = ^o(D(2/, -^-(y, *))>*> + A (D(y, W°(y, t)))y s ; (7.41) 

3 

C 2 (t) = -C 3 (t) = ^((liv ;/ V¥' ; (y. (7.42) 

a 2 (t) = -a 3 (*) = (div y W°(y,t)) Yf , B A = /x (D(y, W°(y)) y/ , (7.43) 
«?(*) = -«2(*) = (div y W e (y,t)) Yf , B{ = ^<B(y, W e (y)) Yf , (7.44) 

B e 2 (t)=^(B(y,^(y,t))) Yf + \ (B(y,W e (y,t))) Ys , (7.45) 

3 

B e = x I + ^{xo^VeV, + x 0s (VeVJ ® e,-, (7.46) 

i=l 

where x = mx / + (1 — m)x 0s . □ 

Lemma 7.4. Tensors A 2; A3 and A4, matrices B4, B^, Bq, B\, B e 2 , C2 and 
C3 and scalars a 2 , 03, a? and a 2 are well-defined and infinitely smooth in 
time. The matrix B® is strictly positively defined. 

If a porous space is connected, then the symmetric tensor A2 is strictly 
positively defined. For the case of disconnected porous space (isolated pores) 
A 2 = and the tensor A 2 becomes strictly positively defined. 

The proof of this lemma, except statement about a matrix B®, one can 
find in [TO]. All properties of the matrix B® are well-known (see |14j . [8]). 
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